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Abstract 

The electronic structures of four Laves phase iron compounds (e.g. YFe 2 , ZrFe 2 , LuFe 2 and 
HfFe 2 ) have been calculated by the state-of-the-art full potential electronic structure code. The 
magnetic moments collapse under hydrostatic pressure. This feature is found to be universal in 
these materials. Its electronic origin is provided by the sharp peaks in the density of states near 
the Fermi level. It is shown that a first order quantum phase transition can be expected under 
pressure in Y(Zr, or Lu)Fe 2 , while a second order one in HfFe 2 . The bonding characteristics are 
discussed to elucidate the equilibrium lattice constant variation. The large spontaneous volume 
magnetostriction gives one of the most important character of these compounds. Invar anomalies 
in these compounds can be partly explained by the current work when the fast continuous magnetic 
moment decrease at the decrease of the lattice constant was properly considered. This work may 
remind the experimentalists of these “old” compounds and exploration of the quantum properties 
under high pressures are greatly encouraged. 
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I. INTRODUCTION 


Magneto-structural interaction is a fundamental phenomenon in magnetic functional ma¬ 
terials. The martensitic transformation caused by magnetization in magnetic shape memory 
is a direct magnification of this effect.- Because of the huge variation of the dimensions, the 
compounds can be used in sensors and actuators. Meanwhile, this effect can also compensate 
the variation of the volume due to temperatures, as shown in the so-called Invar effect^. It 
was proposed to be connected with magnetism because the plateau of the volume expansion 
starts apparently almost at the Curie temperature below which ordered magnetic moments 
are established. The nearly zero volume expansion is because that the Invar alloy has a 
spontaneous volume magnetostriction large enough to compensate normal thermal expan¬ 
sions due to lattice vibrations as suggested in the so-called magnetostriction model of Invar-, 
one of the many models to explain this phenomenon. 

The cubic Laves phase iron compounds show the Invar effect in stoichiometry which may 
serves as a simpler model system for developing theory of Invar avoiding disorder complexing. 
It excludes the noncollincar magnetism and antiferromagnetic moment as the mechanism of 
Invar anomaly^. ZrX 2 can be used as a hydrogen storage materials 4 because of its suitable 
binding energy at the interstices. It was found that H-induced lattice expansion will cause 
an increase in the magnetic moment. They were found interesting even half century ago. As 
suggested in the early work of Klein et al.— the electron-phonon coupling in the hypothetical 
paramagnetic ZrFe 2 is substantial and the superconductivity transition temperature can be 
as high as 9 K. 

Early work in AFe 2 found that A and Fe atonies are in ferrimagnetic order. Asano^ stud¬ 
ied the phase stability by comparing the total energies of different phases (nonmagnetic, 
ferromagnetic, and antiferromagnetic states of C14 or C15 Laves phases). He concludes that 
Y, Zr and Hf compounds are ferromagnetic C15 Laves phase at the ground state, which 
is in agreement with the experiments. Related properties of Laves phase iron compounds 
(e.g. YFe 2 , ZrFe 2 , LuFe 2 and HfFe 2 ) were studied in the past. Yamada' has calculated the 
high held susceptibility Xhf of ZrFe 2 being 5.8 x 10 -4 emu/mol, and YFe 2 being 5.57 x 1CU 4 
emu/mol, which agree with the experimental values 6.1 x 10 -4 emu/mol and 1.55 x 1CU 4 
emu/mol reasonably well. Wortmann^ showed that the hyperfine held decreases to zero 
at about 40 GPa in LuFe 2 , and 50 GPa in YFe 2 at room temperature by nuclear forward 
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TABLE I: The experimental (exp.) and calculated values by LDA and GGA of the lattice constant 
(ao), spontaneous volume magnetostriction (u s ) and magnetic moment (M s ) of AFe 2 (A = Y, Zr, 
Hf and Lu) compounds. M s was obtained at T= 4.2 K. 


AFe 2 


a(A) 


c 0, 

,(xl0 3 ) 


M s (^b) 

exp. 

LDA GGA 

exp. 

LDA GGA 

exp. 

LDA GGA 

YFe 2 

7.363 

7.04 

7.28 

small 

49 

55 

2.90 

2.57 

3.21 

ZrFe 2 

7.06 

6.84 

7.04 

10 

15 

54 

3.14 

2.38 

3.10 

HfFe 2 

7.02 

6.82 

7.00 

8 

35 

48 

3.36 

2.86 

3.26 

LuFe 2 7.217 

6.93 

7.15 

r-u 

25 

47 

2.97 

2.53 

3.06 


scattering. At low temperature, the loss of magnetism took place at about 90 GPa. Direct 
information of the magnetic moment under pressures was reported by Armitage^. The mea¬ 
sured ^f-’s are —8.2 ± 0.4 and —6.3 ± 0.3 x 10 -4 kbar -1 for YFe 2 and ZrFe 2 , respectively. 
The experiments reported lattice constants, saturation magnetization and Curie tempera¬ 
tures are listed in Table HI 

Density functional theory (DFT) was a powerful tool to explain and predict the magnetic 
moment under pressure. For example, the HS-LS transition of transition metal monoxides 
(e.g. FeO, MnO, etc.) under the hydrostatic pressure as high as about 200 GPa were 
predicted by CohenM. Magnetic transition in these highly correlated insulators is the results 
of competition among the kinetic energy, exchange energy and Coulombic repulsion? 14 '. The 
magnetic collapse in metals on the other hand can be qualitatively understood with the help 
of the Stoner model: In a simplified version of this model, the magnetic state is stable if 
IN(Ep) > 1, where / is the Stoner parameter, which is weakly dependent on the atomic 
distance, while the density of states at the Fermi level N(E F ) decreases as the band width 
increases under the pressure. At a certain critical pressure, the criterion is no longer satisfied, 
then the ferromagnetism cannot be sustainable. 

In this work, we found that there may exist a first order quantum phase transition 
under pressure in these compound, which is similar to the well studied case of ZrZn 2 , MnSi, 
etc. Large volume magnetostriction, which is an Invar character of these compounds was 
presented. We further suggested that in order to show Invar effect, the magnetic moment 
decrease with the volume should scale with that with the temperature. 


3 









II. CALCULATION DETAILS 


The C15 structure Laves phase (space group Fd3m) has two formula units per face 
centered cubic unit cell. The full-potential local orbital minimum basis band structure code 
(FPLO) 12 was used in our calculation. Both the local spin density approximation (LSDA)^ 
and general gradient approximation (GGA) 14 of the exchange correlations functionals were 
used here and the resuts were compared when necessary. The number of /c-points in the full 
Brillouin zone (BZ) is 30 x 30 x 30, which can guarantee the convergence of the total energy 
to microHartree. The scalar relativistic treatment was used where all the relativistic effects 
were included except the spin-orbital coupling. The fixed spin moment (FSM) calculations 
were used to investigate the possible multiple local energy minima with respect to the 
magnetic moment. The spontanous volume-magnetostriction u s is defined in terms of the 
ratio of the equilibrium volumes in the ferromagnetic FM (' Vfm ) and the paramagnetic PM 
state (V pm) 


V, 


w* = 


FM 


V } 


PM 


VpM 


( 1 ) 


III. THE GROUND STATE PROPERTIES 

The calculated lattice constant (a), the total magnetic moment at equilibrium are listed 
in Table E The agreement with the experiments is reasonably good. The lattice constant 
from LSDA is lower than the experimental ones, which is notorious. The GGA results show 
a much better agreement. However, the volume magnetostriction is largely overestimated 
by GGA than LSDA. The systematic tendency is that the lattice constant and magnetic 
moments of Y and Lu compounds are smaller than these of Zr and Hf compounds, where 
one more d electron is added in the latter compounds. 

The DOS at the equilibrium lattice constant of the selected compounds are shown in 
Figure CD Common features of these DOS are quite obvious. The total DOS are mainly 
contributed by Fe’s, the other elements only show noticeable contributions at energy about 
0.5 eV above the Fermi level. The bandwidth of the 3d states from Fe is about 5 eV, which 
is a typical value of itinerant system. These states are responsible for the properties of these 
compounds. Above the Fermi energy, there is relatively high DOSs in the minority spin 
channel. This feature is dominated by the anti-bonding 3d states of Fe. The bonding and 
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FIG. 1: (Color online) Density of states (DOS) and the atomic projection (shaded area) of 
AFe 2 (A=Y, Zr, Hf and Lu) at the ground state. 

antibonding states of the minority are separated by a deep and wide valley with width about 
1 eV near the Fermi level. 



FIG. 2: The band structure and the DOS of the minority spin channel of ZrFe 2 at the theoretical 
equilibrium lattice constant. The rectangular block with height of 0.06 eV highlights the energy 
window around the pronounced DOS. 

The van Hove singularity just above the Fermi level in the spin down channel is closely 


5 
























































































related to the magnetic properties which we are going to discuss, so we explore the origin of 
it here. Because of the similarity of the electronic structure in these four compounds only the 
bands and the DOS of ZrFe 2 are shown here. The zoomed-in band structure and the DOS of 
the minority states of ZrFe 2 are shown in Figure[21 The peak in the minority states just above 
the Fermi level comes from almost non-dispersive states. This energy range is highlighted 
by the rectangular block with height of 0.06 eV. The high DOS around the peak stems from 
the narrow bands in the directions along W — K and W — U. This small dispersive states are 
stemmed from the specific atom arrangement of the compounds. The Fe atoms in the C15 
Laves phase form the so-called pyrochlore structure with corner shared tetrahedrons. If we 
look at the net in the {111} layer, it is the Kagorne net with alternatingly connected triangles 
and hexagons. As already shown by Johnston and Hoffmann? 15 ', the high peaks in the DOS in 
a Kagorne net of iron atoms come from narrow bands with d — tt character. One band tight 
binding calculation by Isodar^ discovers two non-dispersive degenerated states along the 
X — W line. Furthermore, there are two additional non-dispersive degenerated antibonding 
states along all high symmetry directions. These results indicate that the spiking DOS is 
originated from the special geometrical arrangement of the Fe’s. 



FIG. 3: The “fat” band and PDOS of the minority d-state of Fe in ZrFe 2 . The T 2ff states are from 
states with np=-2, -1 and 1, while the E 9 states are from np=0 and 2. 

If we plot the up-resolved “fat” band (wighted contributions of atomic orbitals) and the 
partial DOS (PDOS) of Fe as shown in Figure El it is clear that at T-point the states with 
rip = 0 (d z 2 ) and 2 (d x 2 _ y 2 ) are degenerate and rip = —2 (d XJ/ ), -1 (d^), and 1 (d^) are also 
degenerate. These two groups are denoted by E g , and T 2ff , respectively. It can be shown that 
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the PDOS is divided into two catalogs: One is from the three d orbitals (rri/ = —2, -1, and 
1) and the other is from the rest two d orbitals (ny=0, 2). It is shown that the strongest van 
Hove singularity in the PDOS just above the Fermi level of the minority spin state is solely 
from the T 2g states. It is quite understandable because the T 2g states form quite strong 7r 
binding states along each atomic chain with nearest neighbor interactions. The antibonding 
state of these n orbitals is the source of the quite spiking feature of the DOS in the cubic 
Laves phase compounds as analyzed by Johnston and Hoffman!- 5 -. From another point of 
view, if we rotate the coordinates and put the z-direction to the diagonal direction of the 
cube, then the lattice can be described by a trigonal lattice. There is a one dimensional 
irreducible representation A lfl with basis of d 2 /2 for the point group of Fe site, of which is 
D 3 d- It turns out that the d 2 /2 orbital comes solely from a linear combination of the original 
T 2g orbitals. As discussed by IsodaJ^ by single orbital tight binding calculations, which is 
naturally a one dimensional representation, the antibonding orbitals are non-dispersive in all 
high symmetric directions of the BZ. This implies that the one dimensional representation 
with d 2 /2 orbitals as its basis should give also quite spiking feature in the DOS. 

IV. THE MAGNETIC MOMENT VARIATIONS UNDER PRESSURE AND ITS 
ELECTRONIC CHARACTERS 

Because of the differences of the A atoms, we can naturally expect some differences 
among these compounds. Firstly the lattice constants of these materials are more or less 
determined by the atomic volume of A. Taking the atomic volume, defined by (atomic 
weight/mass density), of the elements: Y= 19.89, Zr= 14.06, Lu= 17.78 and Hf= 13.41 
(cm 3 /mol), respectively, we can see that the lattice constants in Table Hj follow the same 
tendency. 

Secondly their magnetic moments have different behaviors under pressure. The depen¬ 
dence of the magnetic moment on the lattice constants are shown in Figure |4j The corre¬ 
sponding hydrostatic pressures are shown on the upper abscissas. 

Obviously all of them show a decrease of the magnetic moment with the decrease of 
the lattice constant as expected from the itinerant electron magnetism, but the Hf and Zr 
compounds show a more rapid decrease of the moment at a lattice constant around 6.8 A 
(in the vicinity of the equilibrium lattice constant), while the other two show a gradual 
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FIG. 4: (Color online) magnetic moment variations of AFe 2 (A = Y(a), Lu(b), Hf(c) and Zr(d)) 
at different lattice constant. The corresponding hydrostatic pressures at the lattice constants are 
shown at the right of each subfigures. 

decrease at this low pressure. At high pressure, all four compounds show at least one first 
order transition to a lower or zero spin state. The differences are quite understandable by 
examining the differences of the electron numbers of the compounds under the assumption 
that the electronic structure is not so much influenced by the difference of the A atoms. YFe 2 , 
ZrFe 2 , HfFe 2 , and LuFe 2 show basically similar DOS as discussed before. The difference of 
the electron number shifts the Fermi level in these systems. Zr(4d 2 ) and Hf(5d 2 ) have one 
more d-electron than Y(4d 1 ) and Lu(5d 1 ), so the Fermi levels of the former are shifted 
towards higher energy, closer to the pronounced peak of the minority spin DOS as shown in 
Figure Q3 This accounts for the low pressure instability of the moment. 

The four compounds show multi-step magnetic transitions. This process can be under¬ 
stood by the particular DOS of these compounds. Taking ZrFe 2 as an example, the DOS at 
different lattice constants are shown in Figure [5] (a~d). The lattice constants of each figure 
are indicated by the arrows in Figure BJd) with the corresponding labels of (a), (b), (c), and 
(d). 

From Figure 03(a), it is obvious that at the experimental lattice constant the DOS of the 
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FIG. 5: The total DOS and partial DOS of Zi'Fe 2 at different lattice constants. From (a) to (d), 
the lattice constants are 7.08, 6.85, 6.80, and 6.70 A, respectively. In (a), the partial DOS of Fe 
and Zr are also shown. The Fermi level is indicated by the dashed vertical line at E=0 eV. 


up spin, contributed mainly from Fe, has a gradual increase below the Fermi level, while 
the DOS of the down spin has a wide (~ 0.8 eV) dip below and a sharp increase just above 
Ep. Applying pressure will broaden the band and reduce the width of the dip and decrease 
the magnetic moment. Thus the exchange splitting is reduced. The DOS of the up spin 
and down spin moves towards each other. This gradual decrease of the magnetic moment is 
shown in Figure 0(d) between the arrows (a) and (b). The gradual decrease of the magnetic 
moment continues until the Fermi level passes through the high DOS peak of the minority 
spins, seeing Figure [5jb) and (c). Then the magnetic moment is rapidly reduced, as shown 
in Figure IU(d) when the lattice constant is between 6.85 (arrow (b)) and 6.80 (arrow (c)) A. 

Very small magnetic moment is developed under the high pressures in all of these com¬ 
pounds, as can be observed in Fig. 0J This different behaviors when approaching the 
quantum phase transition can be of interest both experimentally and theoretically. For the 
second order phase transition, the effect of fluctuation was shown to lead to novel electronic 
ground states in magnetic metals such as magnetically mediated superconductivity, partial 


9 

























or quadrupolar order and non-Fermi liquid phases. For the first order phase transition, as 
summarized by Pfleiderer 1 -, it is interesting for a number of reasons: It can drive novel elec¬ 
tronic states, novel types of low lying excitations, or signal the existence of subtle quantum 
correlation effects. In general, peaks of DOS (van Hove singularities) near the Fermi level in 
all real materials result in a ragged free energy landscape. The shape of the DOS thus has 
a connection with the order of the quantum phase transition. It is necessary to give some 
hints about the order of the transition based on our calculations of electronic structures. By 
the simplified Stoner model, the magnetic free energy of the system in the rigid band model 
is expressed by— 

1 

E{m) — / A ^(m')dm' - Im 2 , (2) 

Jo 4 

where A£(m) is exchange splitting as a function of magnetic moments m, and I is the Stoner 
parameter. The A £ (m) can be expanded as a power series of m: 



A£(m) = aim + a 3 m 3 + a 5 m 5 • • • 

( 3 ) 

where 


( 4 ) 


a 3 = i(3A^iVr 5 - iy 3 Af 4 ) 

( 5 ) 


a 5 = 4(105iV 2 4 iV 1 - 9 - lOSlVglVflVr 8 

5! 

TlOlVflVf 7 + 1 5N 4 N 2 N{ 7 - N 5 N{ 6 ). 

(6) 

Ni is defined as the (i 

— l)-th order derivative of the density of states at the Fermi level with 

respect to the energy 

Then the free energy is 



E ( m ) = i(«i - tt)™ 2 + 7«3ra 4 + J:a 5 m 6 

2 2 4 o 

(7) 


The stability of the phase can be discussed in line with Landau’s theory of second order phase 
transitions. Magnetic instability is necessarily given by the condition that a[ = a\ — J < 0, 
which is equivalent to the Stoner criterion IN(Ep) > 1 by considering Equ. ()4f) . 

The necessary condition to have a first order transition is cq — / > 0, a 3 < 0, and a*, > 0 
if higher order terms than m 5 are neglected in Equation ([7j). This means the DOS at the 
Fermi level should be sufficiently small (the Stoner criterion is not fully satisfied) and the 
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curvature of the DOS at Ep is positive and large, so that N s is positive and large enough 
to give negative 03 , otherwise, if IV 3 < 0, a 3 is definitely positive. These first two conditions 
require that the Fermi level is at a narrow valley of the DOS. 

Direct FSM calculation results and the corresponding DOS to analyze the transition were 
added to the above qualitative analysis. The first example is ZrFe 2 , which shows the first 
order transition to the non-magnetic state. The FSM energy curves are shown in Figure [ 6 ] 
at lattice constants around the transition point. The E(m ) curves at a=6.30 A are enlarged 



FIG. 6: The FSM curves of ZrFe 2 at the lattice constants around 6.30 A. The inset shows the 
enlarged curve at the lattice constant a=6.30 A. It clearly shows that magnetic and nonmagnetic 
solutions coexist at this lattice constant. The data in this figure are obtained with 3107 k-points 
in the IBZ. 

in the inset. It clearly shows two energy minima at m=0 and m=0.085 /ig/Fe. The DOS of 
the related nonmagnetic and magnetic solutions are shown in Figure CD It is clear that the 
Fermi level (the dashed vertical line in the figure) is at a dip (between two peaks marked 
by two ellipses) of the nonmagnetic DOS. At the magnetic solution, the two subbands are 
shifted against each other as shown by the dashed horizontal arrows. 

The other example is YFe 2 where the magnetic transition is the second order. The FSM 
curves are shown in Figure 0 The energy minimum moves to zero when compressing the 
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FIG. 7: The DOS of nonmagnetic state (dashed lines) and ferromagnetic state (red lines) of ZrFe 2 
at a=6.30 A. The horizontal dashed arrows show the relative shift of the DOS of the up and down 
spin subbands. The two ellipses indicate the two peaks around the Fermi level which cause the 
first order magnetic transition. 

lattice as shown in the figure. The energy curve at a=5.99 A is zoomed in and shown in 
the inset. The FSM energy difference of small magnetic moments reaches the accuracy limit 
guaranteed by the code. This is the reason that we should use the DOS in order to discuss 
the possible magnetic solutions. The DOS of nonmagnetic and ferromagnetic states are 
shown in Figure [9] It is clear that the “valley” character around the Fermi level is missing 
compared with Figure [3 Rather, Ep is situated at a plateau which cannot have more than 
one magnetic solutions. The other two compounds, LuFe 2 and HfFe 2 show similar second 
order transition. 

V. THE MAGNETOSTRICTION INVAR MODEL 

Understanding the Invar effect, however, has been a problem for half of the century. More 
than twenty different models have been published in the past 50 years for explanation of the 
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FIG. 8: The FSM energy of YFe 2 at lattice constants around a=6.00 A. The inset shows the 
zoomed-in curve at the lattice constant of 5.99 A with an error bar of 0.01 pHartree. 

Invar effect. A general review about the Invar effect can be found, for example, in handbooks 
edited by Buschow and Wohlfarth- 2 , and references therein. One model called 27 -modekhi is 
based on the hypothesis of Weiss that there exist two separated energy minima with different 
volumes and magnetic states: High spin (HS) at high volume and low spin (LS) at low 
volume states. First principle calculations of Fe 3 Ni by EntePl and other authors supported 
the 27 -model. Entel argued that the special position of the Fermi level in the minority band, 
being at the crossover between nonbonding and antibonding states, is responsible for the 
tendency of most Invar systems to undergo a martensitic phase transition. Two minima 
binding curves should lead to some discontinuity (a first order transition) in the pressure 
dependence of certain physical properties, such as volume, magnetic moment etc., but this 
kind of discontinuity has never been observed in Invar alloys as far as we know. This gives 
an obstacle in applying the 27 -model to explain the Invar effect. 

The HS-LS transition can also be continuous and it is in the Invar alloy like ZrFe 2 and 
HfFe 2 as in Figured This point can be clearly illustrated by our FSM calculations. In the 
FSM energy curves, the energy minimum shifts to the lower magnetic moments as the lattice 
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FIG. 9: The DOS of the nonmagnetic state (dashed lines) and the ferromagnetic state (solid lines) 
of YFe 2 at a=5.99 A. The horizontal dashed arrows show the relative shift of the DOS of the up 
and down spin subbands. The Fermi level is shown by the vertical dash. 

constant is decreased as in Figure [TUI Here the FSM energy curves of ZrFe 2 is taken as an 
example. The quite flat FSM energy curves, which means a large spin susceptibility, near 
the transition region, because the average DOS at the Fermi level is large. The reciprocal 
susceptibility, x~m = E" (M), is given by^ 

Xm = hs (2^e// — -0; (8) 

where I is the Stoner parameter. Thermal excitations cause loss of the magnetic moment 
leading to a magnetic transition from the HS state to the LS state. Therefore, increase of 
the temperature leads to gradual loss of the spontaneous volume expansion associated with 
the ferromagnetic state. This gradual process, contrary to the two states (HS and LS) in 
some Invar alloy (e.g. FeaNi), will not cause any discontinuity in the pressure dependence of 
physical properties. In our compounds ZrFe 2 and HfFe 2 the gradual decrease of the magnetic 
moment is the essential difference, compared with the discontinuity present in a typical Invar 
system as FeaNi. 

As an important character of Invar alloy, the spontaneous volume magnetostriction is 
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Magnetic moment(|a B /unit cell) 

FIG. 10: The FSM energy curves of ZrFe 2 near the HS-LS transition regions. 

calculated by Equ. flTI) . The results are listed in Table [U together with the experimental 
data available^. The theoretical values agree with the experimental ones in the sense that 
they are at the same order. The overshooting of the spontaneous volume magnetostriction 
(ca s ) can partly be from the non-vanishing local magnetic moment above the transition 
temperature in the experiments, while in our model it is in a Pauli paramagnetic state where 
the spin moment is zero. The cure for this problem requires a more realistic treatment of the 
paramagnetic phase. It has been shown that a noncollinear— or a disordered local moment 
(DLM) 24 ' 25 model gives a better agreement with the experiments. Nevertheless, the results 
presented here show the major characteristics of Invar alloy: Compared with the compounds 
where no Invar anomaly is observed, the spontaneous volume magnetostriction uj s is larger. 
In typical Invar alloy, such as Ni 35 Fe 6 5 and Fe 7 2 Pt 28 , kh(10~ 3 ) = 18 and 14.422, respectively. 
At the same time, we see that the values of co s of YFe 2 and LuFe 2 are also large. Why do 
they not show Invar anomalies? 

Take a simple ansatz of the temperature dependent volume of a magnetic solid below its 
Curie temperature (T c ), 

V(T) = V 0 + a(T - T C )V 0 + F m (m(r)), (9) 

where Vo is the volume at T c , a is the “non-magnetic” thermal expansion coefficient from 
phonon and electron contributions, and F m (m(r)) is the magnetic contribution to the vol¬ 
ume variation, with normalized magnetic moment m = M/M s and at normalized temper- 
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ature r = T/T c . M s is the saturation magnetization at T=0. The thermal expansion now 
reads, 


dV(T) 

V 0 dT 

dV m (m(r)) 

= a + v v ” 

V 0 dT 

(10) 


M s dV m (m) dm(r) 

= a -\ - 

T c V 0 dM dr 

(11) 


M s dmir) 

= a -\ -. 

" T c V 0 dM/dV m dr 

(12) 


In order to have a zero thermal expansion — ^ } = 0, we require that 


dM M s dm{r) 
dV m T c Voa dr 


(13) 


This shows that the M(V m ) curve should follow the same behavior as m(r), scaled by a 
factor of — T A y oa • As we know that the temperature dependent magnetic phase transition 
is of second order, so m(r) is a continuous function of r. Thus M(V m ) should also be 
continuous. On the other hand, the large decrease of the moment should take place near 
the equilibrium volume at T c because our reference point is T c . In order to show the Invar 
anomaly, the rapid decrease of the magnetic moment should be near the equilibrium lattice 
constant at ambient conditions. This requirement excludes the Y, Lu compounds from Invar 
alloy, where the decrease of the magnetic moment begins too far away from the equilibrium 
volume. Doping of suitable atoms which shifts the Fermi level to the proper place can make 
this transition meet the requirement. 


VI. CONCLUSIONS 

To conclude, we studied the electronic structure and the magnetic moment behaviors of 
four cubic Laves phase iron compounds. The magnetic moment is found to decrease when 
the lattice constant is decreased, and finally disappears. The way of the magnetic moment 
approaching zero can be continuously and discontinuously depending on the geometrical 
characters of the density of states. It can be understood by the Landau’s expansion of the 
magnetic free energy. Invar anomalies in these compounds can be partly explained by the 
current work when the fast continuous magnetic moment decrease at the decrease of the 
lattice constant was properly considered. 
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